o 
o 



(N 



7j 



C/3 

o 

Ph. 



X 



Chaotic Emission from Electromagnetic Systems Considering 

Self-Interaction 

Fernando Kokubun 

Departamento de Fisica, Universidade Federal do 

Rio Grande, 96201-900 Rio Grande, RS, Brazil 

Vilson T. Zanchin 
Departamento de Fisica, Universidade Federal de Santa Maria, 



^ ■ 97119-900 Santa Maria, RS, Brazil 

^ ! Abstract 



The emission of electromagnetic waves from a system described by the Henon-Heiles potential is 



'^ • studied in this work. The main aim being to analyze the behavior of the system when the damping 

c/3 . term is included explicitly into the equations of motion. Energy losses at the chaotic regime and at 

'^ I the regular regime are compared. The results obtained here are similar to the case of gravitational 

Ph. 

waves emission, as long we consider only the energy loss. The main difference being that in the 

(N 

K* ■ present work the energy emitted is explicitly calculated solving the equation of motion without 

■^ _ 

^^ , further approximations. It is expected that the present analysis may be useful when studying the 

CN . 

^T I analogous problem of dissipation in gravitational systems. 
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I. INTRODUCTION 



Motivation 



The chief motivation of the present work is trying to better understand the effects of 
damping forces (the radiation reaction forces) in radiating systems undergoing chaotic mo- 
tions. Our ultimate interest is in gravitational systems, in particular, in the case of radiating 
chaotic systems. However, due to the difficulties that usually arise during the numerical evo- 
lution of equations from Einstein gravity, we start analyzing the electromagnetic analogous 
problem and shall use the experience acquired here to be able to circumvent those difficulties 
in future work dealing with gravitational systems. 

In a classical field theory, the losses of energy and momenta due to the presence of radia- 
tion reaction forces is of fundamental importance to determine the physical properties of the 
system. Studies on this subject have been done in electromagnetic systems since Maxwell 
has established the foundations of the electromagnetic interaction, and in gravitational sys- 
tems just after Einstein has formulated the theory of general relativity. Even though much 
progress have been done in both cases, there are still some points to be clarified. 

Recent works about emission of gravitational waves from chaotic systems presented sev- 
eral interesting features of such systems. However, some important questions remained 
without answer. Particularly, the infiuence of the damping term to the dynamics of a 
chaotic system is not well understood. A major difficulty in studying the effects of radiation 
reaction in the dynamics of a (chaotic) gravitational system is the necessity of including 
higher order Post-Newtonian (PN) terms. Levin 1,] has shown that at 2.0 PN order, the 



Q. 



two body problem with spin is chaotic, extending previous study of Suzuki and Maeda 
Nevertheless, the effects of a dissipation term become important only with inclusion of 2.5 
PN order. So, in order to describe possible effects of chaotic emission on the detection of 
gravitational waves, it becomes important to consider higher order terms (see for example 

n n 

the comments from Cornish [3] and Hughes |^). 

It has been shown that the amount of energy carried away by gravitational waves in a 
chaotic regime is smaller than in a regular regime [5| (see also J2,|3|). However, this result was 
obtained by brute force method, because in Newtonian gravity, the emission of gravitational 
waves is dynamically unimportant, and in these works the energy emission was considered 



at Newtonian approximation. Thus, knowing the exact manner in which the emission of 
gravitational waves in a chaotic system is affected by the damping term is still an open 
question. The way to find the answer to this question is not as straightforward as we might 
naively think. In Einstein gravity, a damping term appears explicitly into the equations of 
motion for a test particle just after some type of approximation is performed. The exact 
form of the dissipation term depends not only on the coordinate system chosen, but depends 
also on the approximation technique used. This is a consequence of the non-linearity of 
the equations of motion. Moreover, the problem of gravitational radiation reaction usually 
involves enormously complicated calculations and are full of potential sources of errors which 
may lead to results whose physical meaning is difficult to be established. 

On the other hand, the analogous problem of the electromagnetic radiation reaction is 
far easier to be analyzed and quite well understood. Much work on the subject has been 
done since the pioneering papers by Lorentz (3\ and Planck |7|. The relativistic version of 
the radiation reaction force was derived by Abraham js| and lately by Dirac 9], and we 
can say that the effects of radiation reaction force on an accelerated particle, as a classi- 

n 

cal field theory in special relativity, is very well understood (see, e.g., !lCi| and references 
therein). The generalization of Dirac's result to curved spacetimes was done by DeWitt 



and Brehme [jjj, and by Hobbs |l^. When considering the quantum theory, the classical 
electromagnetic radiation reaction force is also soundly based, since it can be obtained by 
taking the appropriate limit of a particular quantum electro dynamical process |l3]- For 
instance, the position of a linearly accelerated charged particle in the Lorentz-Dirac theory 
is reproduced by the h — > limit of the one-photon emission process in QED (See jlJ] and 
Refs. therein). However, the study of chaotic radiating electromagnetic systems found in the 
literature refers mostly to quantum properties of such systems. Its classical counterpart was 
not investigated perhaps because the radiation reaction is really important in microscopic 
systems. 

The similarity between the Abraham-Lorentz theory and the equations appearing in some 
approximation schemes from the general relativistic analogous problem of a radiating gravi- 
tational system (see e. g. Ref. U^]), and the simplicity of the electromagnetic case compared 
to the gravitational case, makes interesting to deepen the study on this subject. Therefore, 
we perform here the analysis of the effects of radiation reaction forces considering a classical 
electromagnetic chaotic system, and in a future work we investigate the gravitational case. 



We expect that the comparison of the results from the present work to future works consid- 
ering gravitational systems, although different in characteristics, shed some light helping to 
better understand the gravitational radiation damping problem, particularly in chaotic sys- 
tems (see e.g. [l6| for the comparison among electromagnetic and gravitational non-chaotic 
damped systems). 

B. The problem 

In order to investigate the effects of radiation reaction on the dynamics of an electro- 
magnetic chaotic system, we consider a charged test particle (it can be a macroscopic test 
particle) of mass m and charge q submitted to an external electrostatic field. In such a case. 



the non-relativistic motion of the test particle is governed by the equation 
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^-r = -Text + Tn — ; — , (1) 

dt ^""^ '^ dt ' ^ ' 

where Fext denotes the external force acting on the charged particle, — ^ is the complete 

dt 

(convective) time derivative of the external force, and Xq = 2q'^ / {3mc^) is the characteristic 
dissipation time, which indicates how efficient is the energy emission. The last term of the 
above equation is the particle self-force which arises due to the emission of electromagnetic 
radiation, and is interpreted as a dissipative force. Accordingly, such a term is usually 
referred to as a damping term, and also as a dissipation term, both of which are used 
throughout this paper. 

Accordingly, the names damping term, or dissipation term are also used. 

The derivation of Eq. (0), and of its relativistic version, with some applications and 
with the interpretation of the dissipation term (and, in particular, of the parameter r) can 
)e found in the classical textbooks 10|, |l9[. In fact, in the original derivation by Lorentz 
a] and Planck [7] (the relativistic version was derived by Abraham pi and Dirac 9]) the 
resulting equation of motion is ma = Fext + mTc^ (da/dt), which leads to runaway solutions. 
A way to avoid such a type of solutions is by replacing the time derivative of the particle's 
acceleration da/dt by the time derivative of the external force, m{da/dt) = dF^xt/dt, into 
this equation, what yields Eq. (Q) as a first approximation to the equation of motion for a 
charged particle. A deeper analysis, however, performed in Ref. [l8| claims that Eq. (^ 
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is the correct equation of motion for a charged particle submitted to an external force Fext 



fsee also Refs. 



20 



|21|). 

A further well known property of Eq. (^ is that, for motions within a time interval At 
such that At ^ t^, the radiative effects on the dynamics of the system will be negligible, and 
the last term in Eq. (^ can be neglected. Thus, in order for the effects of the damping term 
to be noticeable, the time of observation must be large compared to Tq. This is equivalent to 
say that the effects of dissipation will be important only for situations in which the external 
force is applied for a time interval At much larger than the dissipation time Tq, At ^ r^. 



These conditions were both taken into account in our simulations (see Sec. IIIB|) . Hence, the 
system we are analyzing here can be interpreted as the analogous to the case of an orbiting 
test body in a weak gravitational field, but considering explicitly the damping term. 

C. The structure of the paper 

In the following section we write explicitly the equations of motion for a test charged 
particle in the Henon-Heiles potential, by assuming a non relativistic motion. Sec. IIIII is 
dedicated to report the numerical results and to their analysis. A brief analysis on the 
relativistic particle motion is done in Sec. IIVI and finally in Sec. |3we conclude by making 
a few remarks and final comments. 

II. HENON-HEILES SYSTEMS 

A. The model 

_We consider an external force Fgxt derived from a Henon-Heiles electrostatic potential 
and work in a non-relativistic regime where Eq. (0) holds (for a relativistic version of 
Eq. ((H) see |l^ |2i|; see also Sec. lIVjl . The choice of such a potential was due mainly to 
its simplicity allied to its dynamical richness, implying for instance chaotic motions, what 
is of capital importance in our analysis. Other interesting point to be mentioned is that 
a potential of the same type was used in a previous work which analyzed the emission of 
gravitational waves p instead electromagnetic waves, and so the results of the two works 



can be compared. Henon-Heiles systems are described by a potential of the form 

V{x,y)=muj — - — + \xy- — ], (2) 



and have been considered in several contexts 



23[ beyond the original astrophysical scenario. 



This potential is basically a perturbed two-dimensional harmonic oscillator. Therefore, u 
may be identified with the oscillatory frequency which, in the absence of the perturbation 



term, is cj = Jk/m, k being a spring constant, for a mechanical system or to = J\Qq\/ma^, 
where Q and q are respectively the source and the test particle charges, for an electric 
system (in CGS-Gaussian units). Parameter a is the characteristic length of the system. 
The characteristic frequency u defines a characteristic period of motion, T = uj~^. 

Without the damping term, and with the usual choice of units [22], m = uj = a = 1, 
and in our case also Q = q = 1 (see below), the chaoticity of the Henon-Heiles system is 
controlled only by its energy E: the system is bound ii E ^ 1/6, being mostly regular for 
the energy range from to nearly 1/10, and being mostly chaotic for E in the range 1/10 
to 1/6. 

In the presence of the damping term, the dynamics of a charged point particle in the 
potential given by equation is governed by the equations 

,2 r 1 



2 ^ 

X + u X = —2 — xy — T„uj 
a 

2 '^ / 2 2 

y + ujy = [x -y 



X + 2- {xy + yx) 
a 



1 , 
y + 2- {xx - yy) 
a 



(3) 
(4) 



Here, working with electromagnetic field and using Eqs. and (j^, we considered the 
effects of radiation damping, comparing long term energy loss between chaotic and regular 
regimes. The energy loss being considered directly into the equations of motion without 
further approximations. The main results are reported and analyzed in Sec. IIHI 

B. Units and normalized parameters 

We present here a discussion about the physical parameters of the present Henon-Heiles 
electromagnetic system. However, let us stress once more that the present model is to be 
considered a toy model, as a laboratory test for our procedures, and not as a test for the 
electromagnetic theory. 



Eqs. Q and Q have three free parameters characteristic to the system under consider- 
ation: The constant a, the characteristic time Xq, and the frequency u. Then we follow the 
standard procedure and choose a new normalized time parameter r and a new normalized 
time variable t given respectively by the relations r = Tq x tu, and t = t{s) x u, where Xq and 
t{s) carry dimensions, while r and t are dimensionless parameters. The constant a, which 
carries dimensions of length, is used to normalize the variables x and y. The usual choice is 
to measure x and y in units of a, which is equivalent to making a = 1 into the equations of 
motion. 

As far as the effects of dissipation are concerned, the important parameter is the rational- 
ized characteristic time r. The contribution of the radiation reaction force to the dynamics 
of the system is proportional to r (see Eq. P). Therefore, the value to be chosen for r has 
to be as large as possible. On the other hand, as shown below, the time of observation (the 
computation time) has to be much larger than r in order for the effects of dissipation being 
noticeable. 

Considering the motion of charged elementary particles, the largest value for Tq follows 
when the test particle is an electron, in which case one has Xq = Tc ~ 6.3 x 10~^^s. If the 
test particle is a proton then r^ ~ 3.4 x 10^^'^s. For macroscopic systems, however, the ratio 
q/m is not fixed and Tg may assume values several orders of magnitude larger than Tg. For a 
charged test body such that q = ae and m = jSnie one has Tq = [a^ / (3)tc. Take, for instance, 
a test body of mass m = lOg and charge g = 1.0 x lO^esu ~ 3.3 x 10^'^C. Then, it follows 
a = 2.1 X 10l^ /? = 1.1 X 10^8, and t^ = 2.5 x lO^^S. 

Regarding to the third parameter of the model, the characteristic frequency uj, one sees 
that it depends also upon the source of the Henon-Heiles potential, being a typical period 
of the system. For an electromagnetic system it is related to the total charge Q of the 
source by a relation of the form, mu'^ = Qq/a^, m and q being respectively the mass and 
the electric charge of the test particle, and a being the characteristic length of the system 
mentioned above. As usual in Henon-Heiles systems, we are free to fix uj as an inverse time 
unit, u = Uq, in such a way that r = TgUg. Thus, if the source of the potential has mass 
m = (3mg and net electric charge g = ae, we have 

a^ 
r = TgUjg = -j-'Te^g ■ (5) 

Therefore, if we think of a specific orbiting particle (r^ fixed), different values of the dis- 



sipation parameter mean different values for cjq, which gives the corresponding physical 
parameters of the Henon-Heiles potential. 

Now we are ready to fix the parameters and to establish some constraints to the physical 
size of our system. Let us consider a system with typical size L (which is essentially of the 
order of the parameter a mentioned above). Being uo'^^ a typical period, we have that LuOq 
is a typical velocity of the system. Such a velocity has to be smaller than the speed of light 
c, i.e., Lujq < c, so that we have the following upper bound for the system size 

L < au;' = c-^ ~ 2 X 10-^=^^ cm. (6) 

^ T p rp 

As we shall see below, physically interesting values of r for which the numerical results 
can be clearly interpreted lay in the interval r G [10^^'^, 10~"^]. Considering such a range for 
r and taking an electron as the test particle, for which a = /? = 1, we obtain the typical 
size of the system as L ^ 2 x lO^^cm for r = 10^^, and L ^ 2 x lO^^cm for r = 10^^'^. On 
the other hand, if the test body has mass m = lOg and charge q = lO^esu {a = 2.1 x 10^^, 
/? = 1.1 X 10^^), we obtain respectively, L ^ 8 x 10~^cm and L ^ 80cm. If the test particle 
is an electron, the typical size of the system is microscopic which is more difficult to be 
managed. Thus, in order to consider a possible experimental setup, it will be certainly more 
feasible to work with a macroscopic test particle. However, the choice of r in the above 
range describes both the microscopic and the macroscopic systems. 

One more issue on the subject of fixing parameters concerns the physical properties of 
the source in the Henon-Heiles system. Namely, the electric charge Q and the characteristic 
length a. As we have seen above, these are related to the characteristic frequency Uq, and 
once we have normalized units through Eqs. (0) and ©, the ratio Q/a^ is fixed as soon 
as we fix the dissipation parameter r. From the above definitions and choices it is found 
Q ~ 1.7 X 10^/3 L/a (CGS), where we assumed that the characteristic size of the source 
is of the same order of magnitude as the parameter L defined in Eq. ©. Hence, in the 
case of the preceding examples it gives the upper limits Q ~ 3.6 esu and 3.6 x 10~^esu for 
the microscopic orbiting particle, respectively with r = 10^^° and r = 10~^. And for the 
macroscopic orbiting particle we find Q ~ 7.1 x 10^^ esu and 7.1 x 10^° esu, respectively, for 
r = 10^10 and r = 10"^ 



III. NUMERICAL SIMULATIONS AND RESULTS 

A. Methodology 

Equations Q and (j3)) were solved numerically. For the sake of comparison, we initially 
used two different numerical methods: a fourth order Runge-Kutta with fixed stepsize and 
a Runge-Kutta with adaptive stepsize |2,4l]. Also, we used MATHEMATICA built-in proce- 
dures for solving Ordinary Differential Equations. 

At first, we integrated Eqs. and (0} considering no dissipation term, i.e., with r = 0. 
The initial conditions were generated at random, fixing only initial energy and choosing 
X = at the start. Since in this case the system is conservative, the energy is a constant of 
motion, and its value was used to check and compare the numerical results obtained through 
different methods. No significant differences were observed, so we adopted a fourth order 
Runge-Kutta in our simulations. 

B. Numerical simulations 

The main concern of this work is answering the question: How does the value of the 
dissipation time r affect the dynamics of system? In particular, we also want to know 
how much energy is radiated by the accelerated particle undergoing chaotic motions in 
comparison to regular motions. It is expected that a large value of r will strongly affect the 
dynamics, in opposition to small values, for which the dynamics of the system should be 
weakly affected. Nonetheless, it remains to be defined what values of r can be considered 
large and what are small ones. For comparison we made simulations for several different 
values of r and same initial conditions. After a few tries, we have chosen five particular 
cases to analyze in more detail. The chosen values are r = 0, 10~^, 10~^, 10~^, and 10~^°, 
with initial energy E = 0.12 and the same initial conditions for all of the five cases. 

The results can be seen in Fig. Q where we plot the Poincare sections for each value of r. 
Each one of the graphs represents the resulting section for only one orbit, corresponding to 
the particular initial conditions we have chosen. Notice that except for r = 10~^, all other 
Poincare sections look very similar, suggesting that r ~ I0~^ or greater are to be considered 
large values, and r is to considered small if its value is of the order of 10~^ or below. 

In order to better understand the behavior of the Poincare sections we evaluate the 
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FIG. 1: Poincare sections for E = 0.12 and several values of r. All sections with the same initial 
conditions. 

percent amount of radiated energy AE as a function of time in each case shown in Fig. Q] 
The resuhs are seen in Fig. |21 where we plot the graphics of AE x t for each case. Note 
that these graphs are the actual data points, and not fits adjusting the data. For instance, 
the lines appearing in the first four graphs of that figure are the result of plotting the set 
of points obtained numerically for each one of the particular orbits chosen to be analyzed. 
Such straight lines indicate that energy emission rate is constant, and that the total energy 
of the system decreases linearly with time. This is so for small dissipation times r, while 
for higher values of r the energy loss rate is not constant with time (see Fig. ^. It is 



seen that for r = the variations in the energy are exceedingly small 
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-11 



and 



look like random variations. This is surely not an effect of dissipation, because the total 
energy dissipated during the integration time is essentially zero. These random variations 
are caused by numerical inaccuracy, as it can be inferred by comparing this to the other 
cases with t 7^ 0, where the energy variations are much larger and systematic, causing the 
energy to decrease monotonically with time. For instance, for r = 10~^ the total energy 
dissipated during the integration time reaches nearly 20 % of the initial value, so that at time 
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FIG. 2: Percentage of energy loss for five values of r, as indicated in the figure. Initial energy is 
E=0.12, and the same initial conditions were used for all cases, as in Fig. ^ 

t ^ 2000 the energy of the system is about E ^ 0.096. On the other hand, for r = 10~^ the 
energy variation reaches nearly 0.20% of its initial value in the same integration time, and 
the energy is reduced to nearly E ~ 0.1197, meaning it is almost a constant of motion. Also 
for r = 10^^ and 10~^° the energy variations reach ~ 10^^ % and ~ 10~^ %, respectively. 
Even though these energy variations are quite small, they are about eight (for r = 10~^) 
and six (for r = 10"^'') orders of magnitude larger then in the case r = 0, and yet we can 
see they cause the energy of the system to decrease monotonically with time. 

As a further example of a large (higher) r value, we performed simulations with r = 10~^ 
and with energy E = 0.12, and the results are seen in Fig. IHl The large graphics shows the 
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FIG. 3: Energy emission with high r term as a function of time. The small graphics is the Poincare 
section of the motion. The initial energy is E = 0.12 and r = 10~^. 

percent variations of the energy, and the small graph is the Poincare section for this special 
orbit. Due to the large dissipation parameter, the motion is highly damped and, after some 
time, rest is attained. At time t = 2000 nearly 90% of the initial energy has been carried 
away by electromagnetic radiation. Note also that the amount of energy emitted does not 
vary linearly with time, as it happens for smaller values of r. 

For our purposes here, high values of r, r ^ 10~^, say, are not interesting because the 
dynamics of the system is highly affected and the comparison to the case without dissipation 
becomes difficult (if not impossible) to be done. Then, we considered only r = 10^^, r = 
10^^, and r = 10^^° in our full simulations. 

Once fixed the values of r, the next step was solving numerically Eqs. Q and Q with 
several initial conditions, and for the energy values E = 0.12 and E = 0.14. A set of 
500 distinct initial conditions was generated using a (pseudo) random number generator 
|25l |. and the same set was used for every combination of the controlled parameters, r E 
{0, 10"*^, 10"*, 10"^°} and E G {0.12, 0.14}. For each pair of these parameters we 
performed 500 simulation, with the r = case being included only for comparison purposes. 

Although we performed simulations also for E = 0.09 and E = 0.10, the respective results 
were not considered in our analysis. In such cases, the number of chaotic motions (typically 
less than 10 in 500 simulations) in our set of results was too small for a good statistics, and 
so they would not be useful in comparing chaotic to regular regimes, which is the basic aim 
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FIG. 4: Poincare sections for E = 0.12 and four values of r, as indicated in each graphics, corre- 
sponding to 500 distinct initial conditions in each plot. 

of the present work. 

C. Results and analysis 

Using the results of our simulations, we constructed Poincare section for each case, all 
of them being drawn on the surface x = in phase space. These sections were analyzed 
in order to separate between dynamics with chaotic motions from dynamics with regular 
motions. 

The graphics in Fig. |3] are Poincare sections obtained for E = 0.12, without damping 
term (r = 0), and with dissipation term for r = 10~^°, r = 10~^, and r = 10~^, as indicated 
in each plot. Although the overall aspects are the same, a detailed analysis of individual 
sections reveals different aspects as seen in Fig. |3[ where we plot Poincare sections for 
particular orbits corresponding to three different sets of initial conditions for each value of 
r. Here a very important result is that with the inclusion of a small damping term, the 
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FIG. 5: Poincare sections for E = 0.12 and four distinct values of r, but with the same initial 
conditions, are shown at the same horizontal level. The corresponding r values appear on top of 
each column. 

overall aspect of motions are the same, and in particular, the energy is nearly constant, so 
that Poincare section is still a good tool in order to classify the motions as regular or chaotic. 
From this figure it is also seen the dependence of the dynamics on the initial conditions, 
besides the dependence upon r. It is also worth saying that a few particular orbits, out of 
the 500 initially chosen, were neglected since it was not clear from the obtained Poincare 
sections whether they correspond to regular or chaotic motions (see Table H}. 

Now with our set of simulations already separated into two sets, one set with ordered 
motions only and the other set with chaotic ones, we calculated for each initial condition 
(in each set) a best linear fit to the energy variations I\E = At + 5, determining A and 
B by using standard techniques of linear regression. Then, with such a set of values for A 
and B we determined the mean values < A > and < B >, and their respective standard 
deviations, a a and as, for each regime of motions. The results are summarized in Table 
m where we show just the values oi < A > and a a- These are more important than the 
values of < -B > and as, because they furnish the (time) rate of energy carried away by 
electromagnetic waves. We also show in that table (last column), the resulting number of 
regular and chaotic motions for each pair of values of the initial energy and dissipation time. 
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<A> 


O-A 


Regime 


# 


10-*' 


-9.92 X 10-5 


1.33 X 10-5 


Regular 


282 


£^ = 0.12 




-9.24 X 10-5 


3.99 X 10-6 


Chaotic 


202 


10-8 


-1.02 X 10-6 


1.29 X 10-^ 


Regular 


236 






-9.18 X 10-^ 


5.39 X 10-8 


Chaotic 


260 


10-10 


-1.08 X 10-8 


9.09 X 10-1° 


Regular 


163 






-9.09 X 10-9 


6.90 X 10-1° 


Chaotic 


337 




10-6 


-1.05 X 10-"^ 


1.53 X 10-5 


Regular 


146 


£ = 0.14 




-9.06 X 10-5 


4.76 X 10-6 


Chaotic 


350 


10-8 


-1.06 X 10-6 


1.51 X 10-7 


Regular 


141 






-9.05 X 10-^ 


5.28 X 10-8 


Chaotic 


359 


10-10 


-1.12 X 10-8 


7.97 X 10-1° 


Regular 


96 






-9.06 X 10-9 


7.42 X 10-1° 


Chaotic 


403 



TABLE I: Mean values of the energy loss rates < A >, and the standard deviations a a for each 
energy and r values. The last column shows the number of motions in each particular regime. 

As mentioned above, some orbits are missing because they could not be classified as regular 
nor as chaotic ones. 

With the obtained data, we compared the amount of energy radiated in regular regimes 
with respect to chaotic regimes, and calculated the percent ratio t] as follows 

<A>R-<A>c 



V 



X 100, 



(7) 



<A>c 

where the subscript C stands for chaotic and subscript R, for regular. This is shown in 
Table IHl In all cases the average energy radiated in regular regimes is larger than in 
chaotic regimes. These results are compatible with what was obtained when considering 
gravitational waves emission 2,|5| (see also l|, |3|). 
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v{%) 


0.12 


10-*' 


7% 


0.12 


10-8 


11% 


0.12 


10-10 


18% 


0.14 


10-*' 


16% 


0.14 


10-8 


17% 


0.14 


10-10 


24% 



TABLE II: The percent values of energy emitted by the system in a regular regime with respect 
to a chaotic regime. 

IV. RELATIVISTIC MOTION 



We have also investigated the behavior of the electromagnetic Henon-Heiles system in 
relativistic dynamics. In such a case we solved the equations |liJ, llS| 



dF 



di ' d± r^ d.t \ J 



dt 



7^ dv 



(8) 



where p = jniv, 7 = 1/Vl — v'^ , and F is the external force given hj F = — Vf/, U being 
the potential function given by Eq. (j21). The explicit form of Eqs. (jHl), analogous to Eqs. 
(JHj) and (JH), were used in the numerical calculations. 

The numerical results obtained from the relativistic equation (jH)) were essentially the same 
as in the non-relativistic case. This can be understood noticing that for the bound system 
the particle undergoes a non-relativistic motion, as can be verified by the following facts. 
In the Henon-Heiles potential Q, for the test particle to acquire velocities comparable to 
the velocity of light, its initial energy E has to be large. In the rationalized units used here, 
this means E' ~ 1. However, as shown in Ref. |22|, if E is larger than 1/6 the system is not 
bound, and then in the regime where relativistic effects become important the particle is not 
bound by the Henon-Heiles potential. Therefore, the relativistic regime is not important in 
the present analysis. 
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V. FINAL REMARKS 

Our results show that when we consider exphcitly the effects of radiation reaction force, 
the energy emission through electromagnetic waves in the chaotic regime is smaller than in 
the regular regime, as it was in the case of emission of gravitational waves. 

The ratio of energy loss in regular compared to chaotic motions increases with the initial 
energy of the system, and decreases with the dissipation parameter. Since in Henon-Heiles 
systems the chaoticity increases with the energy, this means that the ratio between the 
energy emitted in regular motions and in chaotic motions grows with the chaoticity of the 
system. 

We recall that, in the gravitational waves case, the simulations were performed at PN 
approximation lower than 2.5PN. The result was that the effect of gravitational waves emis- 
sion is negligible to the dynamics of the system. However, being PN lower than 2.5, in those 
simulations the effects of radiation emission were in fact not fully considered. In our analysis 
of the electromagnetic Henon-Heiles system, these effects are fully considered through the 
radiation reaction force. Another important result is related to the mean life-time of source. 
If we make a prediction considering only regular dynamics its mean life-time may be shorter 
than the prediction from chaotic dynamics. However, in the case of dissipation by emission 
of gravitational radiation a more careful analysis has to be done. 

The numerical procedures and analysis performed in this work will certainly be useful in 
our task of studying the gravitational analogous problem, the one about the gravitational ra- 
diation emitted by a particle undergoing chaotic motion, considering explicitly the damping 
term into the equations of motion (work on this subject is in progress). 
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